Allele-specific miRNA-binding analysis identifies candidate target genes for breast cancer risk

Most breast cancer (BC) risk-associated single-nucleotide polymorphisms (raSNPs) identified in genome-wide association studies (GWAS) are believed to cis-regulate the expression of genes. We hypothesise that cis-regulatory variants contributing to disease risk may be affecting microRNA (miRNA) genes and/or miRNA binding. To test this, we adapted two miRNA-binding prediction algorithms-TargetScan and miRanda-to perform allele-specific queries, and integrated differential allelic expression (DAE) and expression quantitative trait loci (eQTL) data, to query 150 genome-wide significant ( P≤5×10-8 ) raSNPs, plus proxies. We found that no raSNP mapped to a miRNA gene, suggesting that altered miRNA targeting is an unlikely mechanism involved in BC risk. Also, 11.5% (6 out of 52) raSNPs located in 3'-untranslated regions of putative miRNA target genes were predicted to alter miRNA::mRNA (messenger RNA) pair binding stability in five candidate target genes. Of these, we propose RNF115, at locus 1q21.1, as a strong novel target gene associated with BC risk, and reinforce the role of miRNA-mediated cis-regulation at locus 19p13.11. We believe that integrating allele-specific querying in miRNA-binding prediction, and data supporting cis-regulation of expression, improves the identification of candidate target genes in BC risk, as well as in other common cancers and complex diseases.Funding Agency Portuguese Foundation for Science and Technology CRESC ALGARVE 2020 European Union (EU) 303745 Maratona da Saude Award DL 57/2016/CP1361/CT0042 SFRH/BPD/99502/2014 CBMR-UID/BIM/04773/2013 POCI-01-0145-FEDER-022184info:eu-repo/semantics/publishedVersio

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

Using GWAS top hits to inform priors in Bayesian fine-mapping association studies

The default causal single‐nucleotide polymorphism (SNP) effect size prior in Bayesian fine‐mapping studies is usually the Normal distribution. This choice is often based on computational convenience, rather than evidence that it is the most suitable prior distribution. The choice of prior is important because previous studies have shown considerable sensitivity of causal SNP Bayes factors to the form of the prior. In some well‐studied diseases there are now considerable numbers of genome‐wide association study (GWAS) top hits along with estimates of the number of yet‐to‐be‐discovered causal SNPs. We show how the effect sizes of the top hits and estimates of the number of yet‐to‐be‐discovered causal SNPs can be used to choose between the Laplace and Normal priors, to estimate the prior parameters and to quantify the uncertainty in this estimation. The methodology can readily be applied to other priors. We show that the top hits available from breast cancer GWAS provide overwhelming support for the Laplace over the Normal prior, which has important consequences for variant prioritisation. This work in this paper enables practitioners to derive more objective priors than are currently being used and could lead to prioritisation of different variants

Proteomic analysis in lupus mice identifies Coronin-1A as a potential biomarker for lupus nephritis

Background Approximately 50% of systemic lupus erythematosus (SLE) patients develop nephritis, which is among the most severe and frequent complications of the disease and a leading cause of morbidity and mortality. Despite intensive research, there are still no reliable lupus nephritis (LN) markers in clinical use that can assess renal damage and activity with a high sensitivity and specificity. To this end, the aim of this study was to identify new clinically relevant tissue-specific protein biomarkers and possible underlying molecular mechanisms associated with renal involvement in SLE, using mass spectrometry (MS)-based proteomics. Methods Kidneys were harvested from female triple congenic B6.NZMsle1/sle2/sle3 lupus mice model, and the respective sex- and age-matched C57BL/6 control mice at 12, 24 and 36 weeks of age, representing pre-symptomatic, established and end-stage LN, respectively. Proteins were extracted from kidneys, purified, reduced, alkylated and digested by trypsin. Purified peptides were separated by liquid chromatography and analysed by high-resolution MS. Data were processed by the Progenesis QIp software, and functional annotation analysis was performed using DAVID bioinformatics resources. Immunofluorescence and multiple reaction monitoring (MRM) MS methods were used to confirm prospective biomarkers in SLE mouse strains as well as human serum samples. Results Proteomic profiling of kidney tissues from SLE and control mice resulted in the identification of more than 3800 unique proteins. Pathway analysis revealed a number of dysregulated molecular pathways that may be mechanistically involved in renal pathology, including phagosome and proximal tubule bicarbonate reclamation pathways. Proteomic analysis supported by human transcriptomic data and pathway analysis revealed Coronin-1A, Ubiquitin-like protein ISG15, and Rho GDP-dissociation inhibitor 2, as potential LN biomarkers. These results were further validated in other SLE mouse strains using MRM-MS. Most importantly, experiments in humans showed that measurement of Coronin-1A in human sera using MRM-MS can segregate LN patients from SLE patients without nephritis with a high sensitivity (100%) and specificity (100%). Conclusions These preliminary findings suggest that serum Coronin-1A may serve as a promising non-invasive biomarker for LN and, upon validation in larger cohorts, may be employed in the future as a screening test for renal disease in SLE patients

Familial Glucocorticoid Receptor Haploinsufficiency by Non-Sense Mediated mRNA Decay, Adrenal Hyperplasia and Apparent Mineralocorticoid Excess

Primary glucocorticoid resistance (OMIM 138040) is a rare hereditary disease that causes a generalized partial insensitivity to glucocorticoid action, due to genetic alterations of the glucocorticoid receptor (GR). Investigation of adrenal incidentalomas led to the discovery of a family (eight affected individuals spanning three generations), prone to cortisol resistance, bilateral adrenal hyperplasia, arterial hypertension and hypokalemia. This phenotype exacerbated over time, cosegregates with the first heterozygous nonsense mutation p.R469[R,X] reported to date for the GR, replacing an arginine (CGA) by a stop (TGA) at amino-acid 469 in the second zinc finger of the DNA-binding domain of the receptor. In vitro, this mutation leads to a truncated 50-kDa GR lacking hormone and DNA binding capacity, devoid of hormone-dependent nuclear translocation and transactivation properties. In the proband's fibroblasts, we provided evidence for the lack of expression of the defective allele in vivo. The absence of detectable mutated GR mRNA was accompanied by a 50% reduction in wild type GR transcript and protein. This reduced GR expression leads to a significantly below-normal induction of glucocorticoid-induced target genes, FKBP5 in fibroblasts. We demonstrated that the molecular mechanisms of glucocorticoid signaling dysfunction involved GR haploinsufficiency due to the selective degradation of the mutated GR transcript through a nonsense-mediated mRNA Decay that was experimentally validated on emetine-treated propositus' fibroblasts. GR haploinsufficiency leads to hypertension due to illicit occupation of renal mineralocorticoid receptor by elevated cortisol rather than to increased mineralocorticoid production reported in primary glucocorticoid resistance. Indeed, apparent mineralocorticoid excess was demonstrated by a decrease in urinary tetrahydrocortisone-tetrahydrocortisol ratio in affected patients, revealing reduced glucocorticoid degradation by renal activity of the 11β-hydroxysteroid dehydrogenase type 2, a GR regulated gene. We propose thus that GR haploinsufficiency compromises glucocorticoid sensitivity and may represent a novel genetic cause of subclinical hypercortisolism, incidentally revealed bilateral adrenal hyperplasia and mineralocorticoid-independent hypertension

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis

Somatic Mutagenesis with a Sleeping Beauty Transposon System Leads to Solid Tumor Formation in Zebrafish

Large-scale sequencing of human cancer genomes and mouse transposon-induced tumors has identified a vast number of genes mutated in different cancers. One of the outstanding challenges in this field is to determine which genes, when mutated, contribute to cellular transformation and tumor progression. To identify new and conserved genes that drive tumorigenesis we have developed a novel cancer model in a distantly related vertebrate species, the zebrafish, Danio rerio. The Sleeping Beauty (SB) T2/Onc transposon system was adapted for somatic mutagenesis in zebrafish. The carp ß-actin promoter was cloned into T2/Onc to create T2/OncZ. Two transgenic zebrafish lines that contain large concatemers of T2/OncZ were isolated by injection of linear DNA into the zebrafish embryo. The T2/OncZ transposons were mobilized throughout the zebrafish genome from the transgene array by injecting SB11 transposase RNA at the 1-cell stage. Alternatively, the T2/OncZ zebrafish were crossed to a transgenic line that constitutively expresses SB11 transposase. T2/OncZ transposon integration sites were cloned by ligation-mediated PCR and sequenced on a Genome Analyzer II. Between 700–6800 unique integration events in individual fish were mapped to the zebrafish genome. The data show that introduction of transposase by transgene expression or RNA injection results in an even distribution of transposon re-integration events across the zebrafish genome. SB11 mRNA injection resulted in neoplasms in 10% of adult fish at ∼10 months of age. T2/OncZ-induced zebrafish tumors contain many mutated genes in common with human and mouse cancer genes. These analyses validate our mutagenesis approach and provide additional support for the involvement of these genes in human cancers. The zebrafish T2/OncZ cancer model will be useful for identifying novel and conserved genetic drivers of human cancers